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Abstract: Heavy vector-like quarks of charge —1/3, B, have been searched for at the LHC 
through the decays B —)• bZ, hh, tW. In models where the B quark also carries charge under 
a new gauge group, new decay channels may dominate. We focus on the case where the 
B is charged under a 17(1)' and describe simple models where the dominant decay mode is 
B —7- bZ' —)• b{bb). With the inclusion of dark matter such models can explain the excess 
of gamma rays from the Galactic center. We develop a search strategy for this decay chain 
and estimate that with integrated luminosity of 300 fb“^ the LHC will have the potential to 
discover both the B and the Z' for B quarks with mass below ~ 1.6 TeV, for a broad range 
of Z' masses. A high-luminosity run can extend this reach to 2 TeV. 
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1 Introduction 


Massive vector-like quarks exist in many extensions of the Standard Model (SM), e.g. extra- 
dimensional models (both warped and flat), little Higgs theories, and composite Higgs models, 
and they are being actively searched for at the LHC. Because these massive states are vector¬ 
like they need not have the same SM quantum numbers as states in the SM, but in many 
instances they do. We focus on that case here. In particular, we consider massive quarks, H, 
that have the same SM charges as the right-handed bottom quark. 

These new particles can be produced through their QCD couplings and are presently 
searched for through the decays B —>■ hhiZhiWt [1-3]; similarly, heavy top partners are 
searched for in decays T —)• ht/Zt/Wh [3-7]. The present bounds on the B mass vary from 
~ 750 GeV if the decay is purely to Zb, to ~ 900 GeV if the decay is purely to hb. The 
bound is ~ 790 GeV in the Goldstone limit where the branching ratios are B(B ^ Zb) : 
B{B ^ Wt) : B{B ^ hb) = 1:2:1. The bounds can be weakened if the B quark decays to 
alternative final states. In this paper, we devise an LHC search strategy appropriate for one 
such exotic decay and estimate its potential sensitivity. 

The B quark can be part of a larger extension of the SM and in particular could be charged 
under additional gauge groups. Here we consider a simple extension where the B quark, 
which mixes with the SM b quark, carries an additional U{1) charge. Such a scenario has a 
simple realisation within the context of “Effective Z' models” [8]. These models introduce, 
in addition to the massive vector-like quark, a new 1/(1)' gauge group and a scalar to break 
it. Although we focus on the case where only the vector Z' is lighter than the B, our collider 
analysis will be effective provided that one or both of the Z' and the scalar (j) are lighter 
than the B. In Section 2 we describe in more detail the particle content, parameter space, 
and phenomenology of this class of models. We demonstrate that it is natural for the new 
decay chain B' —)• bZ' —)• b{bb), shown in Figure 1, to dominate over the modes that are 
currently being searched for. We also outline other interesting final states, involving SM 
bosons, leptons or missing energy, that can occur in some regions of parameter space and 
which are also interesting to search for at the LHC. 

There may be other states charged under the 1/(1)', and if any are stable and electrically 
neutral they can be a dark matter (DM) candidate. The annihilation products of such a DM 
candidate would be rich in b quarks. This presents an intriguing possibility since it is well 
known that the excess of high energy gamma rays seen coming from the proximity of the 
Galactic center [9, 10] can be explained by a 30 - 50 GeV DM particle annihilating to bb, 
or a heavier DM particle annihilating to a pair of resonances, with mass near 50 GeV, that 
decay to bb. Thus, there is a possible connection between an astrophysical signal in gamma 
rays and a collider search in multi-6 final states. We will discuss the phenomenology of the 
model once DM is added, and we will include as one of our collider benchmarks a scenario 
where the Z' has a mass of 50 GeV. 

Having motivated B —?■ {Z'/(j))b, Z'/cj) bb as a search channel for heavy B quarks 
we propose a new search strategy at the LHC, described in detail in Section 3. The final 
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state contains six b quarks but due to the kinematics may not contain six 6-jets. For this 
reason, and to be conservative, we only require three 6-tags in each event. To further suppress 
background we find it beneficial to place a cut on the total hadronic activity in the event, 
Ht = that scales with the B mass being searched for. 

To maximize our sensitivity over a broad range of B and Z' masses we apply three 
approaches to event reconstruction, which use the hardest four, five, and six jets, respectively. 
A given event is subjected to all reconstruction methods for which it qualifies, e.g. if the event 
has six or more hard jets all three methods are applied. Each reconstruction method first tries 
to form Z' candidates, keeping only those pairs of candidates whose masses are within 10% 
of one another. If Z' candidates are found we then attempt to form B candidates by pairing 
Z' candidates with an extra jet, and again keep only those that are within 10% in mass. The 
six-jet analysis reconstructs Z' candidates as dijet pairs, the four-jet analysis reconstructs Z' 
candidates as single jets with sub-structure, using the A^-subjettiness variable [11], and the 
five-jet analysis reconstructs one Z' candidate as a dijet system and the other as a single 
jet with substructure. For signal events the distribution of {Mz',Mb) pairs has a clear 
concentration close to the expected values. The background distribution, coming dominantly 
from ti and QCD multi-jet, has a different shape, allowing separation of signal and background 
over a broad range of masses. 

In Section 4 we present our results, which show that discovery at the 5a level is possible 
for a broad range of Mz', with Mb < 1250 GeV for 30 fb“^, with Mb < 1600 GeV for 300 
fb“^, and with Mb < 2000 GeV for 3000 fb“^. Accurately modelling the QGD background 
is a fraught enterprise. In a full experimental analysis the background needs to be estimated 
from data, and we describe one approach to doing so in Section 4. By relaxing the number 
of 6-jets required for an event to pass the cuts one can determine the expected shape of 
the {Mz', Mb) distribution for background alone. The normalisation of the distribution can 
be estimated by comparing the total number of events with and without the 6-tags, before 
the analysis cuts requiring B and Z' candidates. We show that this approach works well 
when tested out on Monte Garlo data and propose other sidebands that may be available to 
estimate the QGD background from data. 

2 An Effective Z' Model 

In this section we describe a particular effective Z' model [8] and identify parameter space that 
realizes the phenomenology we wish to study. Although we add a relatively modest number 
of new fields beyond those of the SM, several new interactions are allowed and multiple new 
phenomena can arise. We introduce a pair of vector-like quarks, {B,B^), which are charged 
under a new 17(1)' and also charged under the SM in a similar way to the RH bottom quark, 
i.e. B has quantum numbers (3,l,l/3,—1) under (56/(3), S'f7(2), [/(l)y, 17(1)') and has 
(3,1, —1/3,1). Because the new quarks enter as a vector-like pair, there are no issues with 
gauge anomalies. In addition we introduce a new complex scalar that has charge -|-1 under 
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Figure 1: The BB production and decay process that is the primary focus of our analysis. 


the but which is otherwise neutral. We assume that $ gets a vev that breaks the ?7(l)^ 
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leading to a mass for the 11(1)' gauge field, 


Mz' = gw. 


( 2 . 2 ) 


For the collider phenomenology that interests us, this is the minimal model. If there are 
also vector-like fermions (x^X^) that are neutral under the SM but charged under the 1/(1)', 
they can provide a viable DM candidate, as we investigate below. An analogous setup with 
a vector-like top quark, T, in place of B has been considered in Ref. [12]. 

The QCD cross section for BB pair production depends only on the mass of the vector¬ 
like quarks, but the resultant final states for these pair-production events depend upon the 
sizes of the various possible couplings between the SM and the new sector. In Section 2.1 
we consider these couplings and the mixings they induce. In Sections 2.2-2.4, we study the 
decays of B, Z', and (j). We find that the decay chain that we use for our collider studies, 
B —)• bZ' —)• b{bb), can easily dominate, although the analysis we develop is equally effective 
if R —)■ —>■ b(bb) dominates. We discuss DM phenomenology in models that incorporate 
the (X)X^) fields in Section 2.5. 

2.1 Mixing of B, Z', and 0 with Standard Model fields 

2.1.1 Quark mixing 

If the only interactions of the vector-like quarks were their gauge interactions, there would be 
an unbroken Z 2 parity under which the new fermions are odd. However, the gauge symmetries 
of the theory allow a so-called 4>-kawa interaction, X^Bb‘^, which breaks the Z 2 and allows B 
to decay. Including this Lagrangian term, the B and b masses arise from 


C D -msBB'' - X^Bb^ - VbQHb^ + h.c. 


(2.3) 
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More generally, can couple to a linear combination of (B, and but to be consistent 
with flavor constraints we assume that this linear combination is dominated by h^. Alterna¬ 
tively, we could introduce three copies of the heavy vector-like quarks that couple in a flavor 
symmetric fashion to the SM down-type quarks, but with a hierarchy in the masses of the 
heavy quarks, such that the only sizable effective coupling of the Z' is to the b quark. 

Either way, once <1> acquires a vev it induces B — b mixing. This mixing is largest in the 
RH quark sector. The mass-eigenstate RH quark fields are 

b^ = CRb^ - srB^ , B^ = crB^ + SRb^ , (2.4) 


with the mixing angle determined by 

SR = sin Or 


m 

Mb 


Xw 

y/2MB 


(2.5) 


Here Mr = is the physical mass of heavier eigenstate, and we work in the 

approximation that the mass of the bottom quark can be neglected. 

The mixing in the LH quark sector is related to the RH mixing by 


th 


Mb 

Mb 


tR , 


( 2 . 6 ) 


where as above we denote the physical mass of a field / by Mj. One consequence oi b — B 
mixing is that the coupling t/b differs numerically from the SM bottom Yukawa coupling, 


Vb = 




(2.7) 


where v ~ 246 GeV. 


2.1.2 Gauge kinetic mixing 

Another renormalizable interaction allowed by the symmetries of the theory is kinetic mixing 
between the U{1)' gauge field {b^) and the hypercharge gauge field (H^), 

C D . (2.8) 

This operator allows the Z' to decay to SM fields. If this operator is absent at some high 
scale A (for example, this could be the scale at which SU{2)' breaks to U{!)'), it will be 
generated by B and b loops. Taking Mr to be somewhat above the U{iy breaking scale, we 
can approximate the value of k at the scale Mr by ignoring the quark mixing, giving 


K ~ 


gyg 


log 



(2.9) 


Provided A is not too far above Mr, we expect k r-^ 10 ^ — 10 ^ for g' ~ qy- Significantly 
smaller values of k are possible for smaller g' , or if contributions from additional states 
partially cancel contributions from b and B loops. 
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Working to first order in k, we obtain diagonal kinetic terms and mass terms with the 
field redefinitions 


t Czb^ (2.10) 

t A^ i^c^iczb^ -|- SzZ^') (2.11) 

t (^Cz KSy/Sz)Z^ -|- (yHS\YCz Sz)b^, (2.12) 


where s\y and C]y are sine and cosine of the weak-mixing angle, and the mixing angle 0z 
is introduced to remove mass mixing induced by the kinetic mixing. This mass mixing is 
required to be small by precision studies, and for k <C 1 it is guaranteed to be small unless 
Mz and Mz' are very close. Assuming 0^ <C 1 and using the leading-order result 


0z = 


Ml 


Ml - Ml, 




(2.13) 


the couplings of the Z' to SM fermions can be determined from 


CD g 


swuM^, ^ ' 




\ , SWKMI ry ' 

Ml-Ml ^ 


(2-14) 


to first order in k. In Section 2.3 we consider the competition between quark mixing and 
kinetic mixing in determining Z' branching ratios. 


2.1.3 Scalar mixing 

With the addition of the scalar potential is 


+ Xh\H\^ - ml\^^ + . (2.15) 


The mixed quartic term leads to a mass mixing between the Higgs and (p fields, producing 
mass-eigenstate scalars 

h = Chh-Sh(t> ^ = Ch4> + Shh, (2.16) 


where 


sin 20/j 


m2 

<P 


M\y 

99 ' 

-^>1. 


(2.17) 


determines the mixing angle. 

Scalar mixing leads to corrections to the partial widths of the SM Higgs boson of the 
form T —)> c^Tsm, with the exception of the partial width to b quarks, which is also altered 
by the b — B mixing. At tree level we have 


r(/i ^ hb) 


‘ic\{cLChyb + SLSh>^f 
IOtt 


Mh = c\cl 



1 SR th Xv 
\/2 Cl Mb 


2 

CsM^h -D bb) . (2.18) 


In the absence of scalar mixing, the correction factor is 


4-1-2 



Mb) 
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(2.19) 













and the deviation from the SM resnlt is tiny due to the smallness of Mf,. 

If the Z' is light enough, scalar mixing also induces a new decay mode, 


r(/i ^ Z'Z') 


g'‘^s\ -y/l “ 4 q:^/ (l - 4 q:^/ + 12 q;|,) 

327r az' 


Mh , 


( 2 . 20 ) 


where az' = This could lead to many interesting signatures depending on how the 

Z' decays, e.g. /i —)• 46, h —)• invisible (if Z' decays to DM), or 6, —)> 4£ without a Z resonance. 
Furthermore, the Higgs may be produced in B decays (as discussed in Section 2.2), resulting in 
a final state from BB production with as many as 10 6’s. Exotic Higgs decays, e.g. h —)• ZZ', 
can also be induced by kinetic mixing. The effects of scalar and kinetic mixing on Higgs 
decays have been widely studied in the literature, see for example Ref. [13]. 

Beyond its effects on the Higgs particle, scalar mixing also impacts (p decays. In Sec¬ 
tion 2.4 we consider the competition between the <h-kawa interaction and scalar mixing in 
determining p branching ratios. 


2.2 Heavy quark decays 


As discussed above, the \^Bb^ interaction term breaks the Z 2 parity acting on the new 
fermions and allows the B to decay. At tree level, the possible two-body final states are Z'b, 
Zb, W~t, pb, and hb. 

For decays of B into a vector boson v and a fermion /, the relevant interaction term has 
the form 

C D f'y^iaiPL + 0'rPr)B v^, (2-21) 

and the tree-level partial width is 


T{B Vf) = ^ ((1 -Xf - Xyf - 4:XfXy) 


1/2 


|2 , I |2' 

Oil + |«r| 


1 + Xy — 2Xf + XyXf — 2x^ + Xf 

- -z-- 


( 2 . 22 ) 

6{alaR +apfOi) y/xj • 


Here we define Xf = Mj/M^ and Xy = M^/M^. Neglecting corrections induced by kinetic 
mixing, the relevant couplings for B —?■ Z'b, B —)• Zb, and B —?■ W~t are 


Ur — 

■ -g'cLSL 

^ = -g'cRSR 

(2.23) 

af = - 

e CRSL 

af = 0 

(2.24) 

2cwsw 


e SL 

=0, 

(2.25) 


where sl, cr, sr, and cr describe the mixing in the fermion sector, with the left- and right- 
handed mixings related through Equation (2.6). 

For decays of B into a real scalar s and a fermion /, the relevant interaction term has 
the form 

CDf{f3LPL + PRPR)B s, (2.26) 


- 7 - 









and the tree-level partial width is 

Mb 


V{B ^ sf) = 

X 


((1 -Xf - Xsf - AXfXsY^"^ 

'l\^ + I/5-RI^) (1 + 3^/ “ ^s) + 2 {(^bPr + I^rPl) y/XfXg 


327r 
2 


(2.27) 


with Xs = M'^ /M"^. The couplings needed for B ^ (j)b and B ^ hb are 

^ + SLShVb) Pr = ^ {sLCh>^ - CLShVb) 


P’l = ^ (clS/xA - SLChVb) ^ {sLSh^ + CLChVb) • 


Sr 


(2.28) 

(2.29) 


We allow for the possibility of mixing in the scalar sector, with Sh and c/j determined by 
Equation (2.17). 

The comparison between the various B partial widths simplifies if we neglect scalar mixing 
(s/i —)• 0) and work to leading non-vanishing order in (M^/Ms)^. In this approximation we 
find 


T{B - 

> Z'b) 

T{B - 

(f^b) 

r{B Zb) 

T{B - 

hb) 

T{B - 

> Wt) 


(Acfi)^ 

647r 

(Acr)^ 

647r 

{VbSRf 

647r 

{ybSRf 

647r 

{VbSRf 

327r 


(1 - Xz'f (1 + 2xz') Mb 
(1 - x^f Mb 
(1 — xz)^ (1 + 2xz) Mb 
(1 - Xhf Mb 

\J{l-xt- xw f - 4.xtxw 


(2.30) 

(2.31) 

(2.32) 

(2.33) 

(2.34) 


X [(1 - xtf + xw + xt - 2xw)] Mb- 
In the regime where Mb is much larger than all other masses, we have 

T{B Z'b) ~ T{B (j)b) 

T{B Zb) ~ r{B hb) ~ ^r{B Wt) 


r{B Z'b) ~ 


Acr 

VbSR 


T{B Zb), 


(2.35) 

(2.36) 

(2.37) 


consistent with Goldstone equivalence. 

Our collider studies will focus on the decay of B to Z'b. As shown in the left-hand plot 
of Figure 2, this decay can easily dominate over decays into SM states, due to the smallness 
of Ub- In fact, using Eqn. (2.5), the quantity appearing on the vertical axis can be rewritten 

crX crMb 


SRVb 




as 


(2.38) 

















Figure 2: Left: Contours of T{B —^ Z'b)/T{B —^ hb/Zb/W t). We neglect the masses of 
all SM particles, which overestimates the partial widths into SM states. Right: Contours of 
r{B ^ Z'b)/r{B ^ (j)b). 

which goes to in the A —?■ 0 limit. It is not therefore not necessary for A to be 

large for B —?■ Z'b to dominate. Given ample phase space for the decay, B —)• Z'b dominates 
over decays to SM states for small A, unless (<h) is much larger than Mb- 

The remaining competing decay, B (pb, can be forbidden kinematically by raising M^j, 
above Mb- A light Z' is consistent with M^ > Mb because g' can be taken to be small. 
The opposite scenario is also possible: one can have a light cp with Mz' > if the quartic 
coupling A(^ is small. In this case B ^ (pb can be the dominant decay. The right-hand plot 
of of Figure 2 shows how the ratio r(i? —)• Z'b)fT{B —)> (pb) depends on and Mz' when 
both channels are kinematically accessible. 

li B ^ (pb dominates, the results of our collider studies apply essentially unchanged, 
provided that cp decays dominantly to bb {cp decays are studied in section 2.4). If instead both 
B —)• Z'b and B ^ (pb have sizable branching ratios, the analysis we develop below is flexible 
enough to reconstruct both BB —>■ cpcpbb events and BB —)• Z'Z'bb events, even if Mz' and M^ 
are very different. Two invariant mass peaks at distinct values of Mz'/M^ would be found, 
with reduced strength compared to the case with just one dominant channel. Our analysis is 
not designed to reconstruct BB —?■ cpZ'bb events efficiently, unless the Z' and (p happen to be 
close in mass. 

2.3 Z' decays 

At tree level, and neglecting kinetic mixing, the potential two-body channels for Z' decay are 
bb, bB, bB, and BB, some of which might be kinematically forbidden. Kinetic mixing allows 
for decays into other fermions, including leptons, and decays to bosons. If DM is charged 
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under U{iy and is sufficiently light, there will also be invisible decays of the Z' , as discussed 
in Section 2.5. 

For decays of the Z' into fermions / 1 / 2 , the interaction term 


leads to the tree-level partial width 

hi 2 ) = ((1 - - Vif - 


(2.39) 


(2.40) 


plI' + \pr\) (2 - yi - y2 - {yi - y2f) + 6 {pIpr + p*rpl) y/mm 


where Nc is the fermion color multiplicity and yi ^2 = 

Neglecting corrections induced by kinetic mixing, the relevant couplings for Z' —)■ bb, 
Z' —)• hB/Bb, and Z' —)• BB are 



pf = 9 'si 

II 

= 9 'Sr 

Bb 

bB / 

Bb 

bB 

Pl 

= Pl = -9 SLCL 

PR ■ 

- Pr - 


BB / 2 

BB 

/ 2 


Pl =9 Cl 

PR 

= 9 Cr 


(2.41) 

(2.42) 

(2.43) 


Dropping terms involving and we hnd 

r(z' ^ bb) = 

OTT 

/2 2 2 

r{Z' ^ bB) = T{Z' ^ iB) = (1 _ y^)2 (2 + yj,) Mz' 

loTT 

T{Z' —> BB) = ^ Y^l — AyB [(l + Cr) — ys (l — 6c^ -|- c|j)] 


(2.44) 

(2.45) 


Mz' .(2.46) 


Our collider studies will focus on scenarios with Mr > Mz', in which case Z' —)• bb is the 
only allowed decay among those above. 

Kinetic mixing modihes the Z' widths given in (2.44)-(2.46) and opens up new Z’ de¬ 
cay modes. If only kinetic mixing is present, the couplings of Z' to SM fermions can be 
summarized as 



cw 


( F/ Qf [1 <lvyz] \ 

V J 



enQf 

Cw 


f l-cl,yz \ 
V ^-yz ) 


(2.47) 


where we work to leading order in k and where yz = These couplings can be used 

with Equation (2.40) to calculate the Z' partial widths into SM fermions induced by kinetic 
mixing. For fermions that can be approximated as massless, the result simplihes to 


T{Z' ^ //) 



^ cl^yzY + 

Qf [1 c^yz] 

n 

247rc^ 

[ “ 

-yz)^ 

1 


(2.48) 


10 
















Figure 3: In the region below the lower orange line (upper blue line), the total width of the 
Z' induced by kinetic mixing (in the absence of quark mixing) is less than 10% (50%) of the 
partial width into bb induced by quark mixing, given by Eqn. (2.44). 


Kinetic mixing also opens up decays of the Z' to boson pairs, if kinematically allowed, with 
partial widths 


T{Z' WW) 
T{Z' Zh) 


(1 - (1 + 2Qyw + 12y^) 

-o-o- 

1927rc^ (1 - yzf 

—^\/(i - y/i - yz)"^ - ^vhvz 

lyzTTC^ 


(1 - j/fe + yz)'^ + 8yz 

(1 - yzf 


Mz' 


(2.49) 

(2.50) 


If present, scalar mixing modifies T[Z' —)• Zh) and, for sufficiently light (j), induces a partial 
width for Z' ^ Z(j) . 

Large values of n allow for abundant Z' production through its couplings to light quarks. 
The Z' can then decay to leptons, and LHC constraints on dilepton resonances potentially 
become relevant [14]. For smaller k the Z' is mainly produced through its interactions with 
b and B quarks, but interesting leptonic signatures can still be induced by n, e.g. one or two 
dilepton resonances produced in association with 6-jets. 

In Sections 3 and 4 we focus on the case where the dominant decay is Z —)• bb. To 
estimate what values of k are consistent with this scenario, we consider the ratio 


Rz' — 


T{Z') with only kinetic mixing turned on 
r(Z' —7- bb) with only quark mixing turned on’ 


(2.51) 


which depends on Mz' and on the quantity 


5'4 


%/2 


/ kMb \ 
\XsrMz’ ) 


(2.52) 
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If we require Rz' to be small, we get the relatively weak constraints on k shown in Figure 3. 
Taking Mb = 1 TeV, A = 1, and sin^^ = 0.5, Rz' < 0.1 implies k < 9 x 10“^ for Mz' = 
50 GeV and k < 9 x 10“^ for Mz' = 500 GeV. 


2.4 (j) decays 

In our discussion of </> decays we will consider the effects of scalar mixing, but we will neglect 
kinetic mixing. If the scalar mixing vanishes, then at tree level, the potential two-body 
channels for (j) decay are Z'Z', bb, BB, bB, and Bb. Because we are mainly interested in how 
(j) will decay if it happens to be produced in B and B decays, we will take for this 

section, kinematically forbidding decays to BB, bB, and Bb. 

Scalar mixing allows the (j) to acquire the decay channels of the SM Higgs. For any decay 
channel X open to a SM Higgs of mass M^, excluding channels involving b quarks, we have 

T{(f) ^ X) = s\T{h ^ X)\Mh^M^- ( 2 . 53 ) 


The decay width to bb depends on the quark mixing. Working to leading order in M^, the 
tree-level width is 


F((/. ^ bb) 


^cj^jshyb - ChSiX)"^ 

IGtt 


■M^ = 


3 

Stt 






The remaining two-body, tree-level partial widths are 


F((^ —)• hh) 
F(0 ^ Z'Z') 


^Vr^4^(l + 2z,)2(^ + ^' 

327r \w V . 

g'‘^cl yi - 4.ZZ' (1 - 4:Zz' + 122:|,) 

327r zz' 




(2.54) 


(2.55) 

(2.56) 


where Zh = M|/M| and zz' = 

If the heavy B quarks decay mainly to (j)b, the results of Section 4 will apply when (p 
particles decay mainly to bb. Taking Mb = 1 TeV and w = Mfj), we show in Figure 4 the 
parameter space where p ^ bb dominates. When (p is sufficiently heavy to decay to W, Z, 
and h pairs, the mixing in the scalar sector must be very small for bb to dominate over these 
modes. We assume Mz' > M^/2 to make Figure 4, but (p —?■ Z'Z' can easily be the most 
important decay mode if it is kinematically accessible. 


2.5 Dark Matter 

In this work we focus mainly on the LHG phenomenology of the B and the Z'. However, our 
model, over part of the parameter space, also provides a natural explanation for the excess 
of high energy gamma-rays seen coming from the proximity of the Galactic center, the so 
called Galactic Genter Excess (GCE), or Gooperon [9, 10, 15-25]. The spectrum of the excess 
photons is well fit by a 30 - 50 GeV DM particle annihilating directly to bb, as well as by a 
10 GeV DM particle annihilating to r’s. It may also be fit by cascade annihilations of DM 
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M^[GeV] M^[GeV] 

Figure 4: Contours of Br{4> —)> bb) = 0.9 (solid) and 0.5 (dashed), for \sh\ = 10“^ (blue), 
10“^ (red), 10“^ (green), and 10“^ (orange). The left and right plots have Sh > 0 and 
Sh < 0, respectively. We take Mb = 1 TeV and w = M^, and we assume that (j) —)> Z'Z' is 
kinematically forbidden. 


to light mediators which in turn decay to pairs of SM particles [26-29]. In particular, the 
spectrum of the GCE is better fit for annihilations of the form xx {bb){bb) than 

for direct annihilations to 5’s if ~ 30 GeV + Mz'l‘2, [29], e.g.{M^, M^) = (60, 50) GeV. 

We introduce a pair of vector-like fermions, XiX^i with charges Qy^ and —Qy under 
the U{1)' but no SM charge. Provided Qy ^ 0, these fermions are stable at the level of 
renormalizable interactions. Recall that we have normalized the U{iy gauge coupling g' so 
that $, B, and B‘^ have charges -|-1, —1, and -|-1. If Qy is not an integer, an unbroken global, 
abelian symmetry guarantees the stability of X) Even if x, x'^ are not absolutely stable, 
they can easily be stable on cosmological time scales if any non-renormalizable operators that 
induce their decays are generated at the Planck scale or some other very high scale. Provided 
Qy 7 ^ ±1, ±2, there are no operators at dimensions five or six leading to x decays. 

For My > Mz' or My > {Mz' + M^)/2, the annihilation processes xX Z'Z' or 
XX are accessible. This allows for a secluded DM scenario [30], in which the couplings 

that determine the relic abundance are independent of those that determine the DM’s coupling 
to the SM. Focussing on xx Z'Z'^ the non-relativistic DM annihilation rate is 


<yVxx-^Z'Z' 


g'X (, vi.y-'Vi vi.y' 

167rM2 V Ml) \ 2Ml) ' 


(2.57) 


For masses that fit the GGE the correct relic abundance is achieved for g'Qy ~ 0.2. We have 
checked this and other results from this section using micrOMEGAs [31]. If xX is also 

a relevant annihilation channel, slightly smaller values of g'Qy work. 

If My is too light to annihilate into final states involving Z' and (p, the correct relic 
abundance can still be achieved through xX bb, mediated by s-channel Z' exchange. 


- 13 - 
















Neglecting mixing in the LH quark sector, the non-relativistic DM annihilation rate for this 
process is 


3s%g»Ql 

2 ,r (Af|,- 4 Af 2 )H MV V 

With the help of Eqn. (2.5), it is useful to rewrite this as 

3X^Ql Ml / 4M2 \ / Ml \ / Ml \ 

^\x^b-b - ^4 y Ml) V Ml) V 4M2 J • 


(2.58) 


(2.59) 


Unlike the case where the relic abundance is set by XX Z'Z' l(f>Z' , achieving the correct 
relic abundance through xX bb requires that Mb not be too large. The annihilation rate 
is resonantly enhanced for Mz' close to 2M^, but the correct relic abundance can also be 
obtained far off resonance. For example, taking Mb = 1 TeV, = 40 GeV (as preferred 
for the GCE), and Mz' = 250 GeV, we need — 4. Taking A = 1 and maximal mixing 
in the RH quark sector, we get g' = {XMz')/[I^ srMb) = 1/4, and the coupling of the Z' to 
DM is not too large: g'Q^ — 1- 

For M^ < Mz '!the presence of DM coupled to the Z' opens up an invisible decay 
mode with partial width 

r(Z' -V xx) = (1 + 2yx) Mz' , (2.60) 


where = Ml/Ml. The invisible width can easily dominate over the width into bb, 
Eqn. (2.44). In this case BB pair production at the LHG can lead to bb + events targeted 
by standard SUSY searches [32, 33]. 

Because the nucleus has no net 6-charge, direct detection rates are highly suppressed in 
the absence of kinetic mixing. Kinetic mixing leads to a spin-independent coupling of DM to 
the proton, and to a cross section per nucleon 

2 /ioi\2 \ 2 , .o /rr.^_'ir\4 

(2.61) 


cr 


131 


g'Qx 

0.2 


K 

10-4 


50 GeVy 
Mz' ) 


X 10 46 cm^ 


where we normalise to scattering off xenon. For My. ~ 50 GeV, LUX has probed down 
to cr ~ 8 X 10-4® cm^ [34]. Parameters chosen to explain the GCE in the secluded DM 
scenario thus require k < 3 x 10-4 to evade direct detection. Values of k this small are not 
unreasonable, especially given that g' can be small. Taking n to be given by Equation (2.9) 
with the log set to one, the constraint is satisfied for g' ~ 1/20, which requires Qy ~ 4 for 
the relic abundance. The xX bb explanation of the GCE is consistent with values of Mz' 
larger than those preferred by the secluded DM explanation, meaning that LUX constraints 
can be satisfied with larger values of k. 

We have been assuming that x and x^ form a Dirac fermion of mass My , but it is possible 
that the mass eigenstates are Majorana fermions. For example, if Qy = —1/2, the interactions 


£ D + h.c. 


(2.62) 


14 












are allowed, leading to Majorana masses when <1> gets a vev. If these Majorana masses are 
much smaller than the Dirac mass, the relic density calculation does not change much, but the 
cross section for direct detection is dramatically reduced. Larger values of and/or A^c can 
change the phenomenology in various ways, e.g. scalar mixing can induce a Higgs-mediated 
contribution to the cross section for direct detection, final states involving 4> can become more 
important for the relic abundance calculation, and (p can potentially decay invisibly to DM. 

3 Searching at the LHC 

Traditional searches for heavy vector-like B quarks have focused on decays to SM bosons and 
quarks [2, 35, 36]. As we have seen, the presence of Z' and p (and x if DAI is included), can 
significantly alter the phenomenology. Which of the many possible search channels dominates 
depends upon the masses of the new particles and upon the relative sizes of the various 
mixings, namely kinetic mixing, quark mixing, and scalar mixing. 

We will consider the situation where the dominant decays are B —)• Z'b followed by 
Z' —)• bb. As discussed in Section 2, B —)• Z'b tends to dominate for > Mb > Mz', 
unless (<h) is much larger than Mb (see Figure 2), while Z' —?■ bb dominates for Mb > Mz' 
and sufficiently small kinetic mixing (see Figure 3). It will be possible to infer from our final 
results the effect of branching ratios smaller than one. If B decays to both Z'b and pb our 
analysis would find both resonances but at reduced significance, as long as both Z' and p 
decay to bb. 

The sizeable QCD production rate of BB, shown in Figure 5, makes our primary channel 
of interest pp —)■ BB —)> {bZ'){bZ') —)> b{bb)b{bb), which is not presently being searched for. 
Before describing in detail the search strategy we advocate, we briefly discuss other interesting 
channels that are worthy of investigation. 

Although their couplings are suppressed by the quark mixing angle, the Z' and p can 
be singly produced in association with b quarks, which may be forward boosted. If these 
states decay to bb, their existence is probed by LHC searches for bb resonances produced in 
association with b quarks [37]. 

With kinetic mixing the Z' will have a di-leptonic branching ratio, but unless k is suffi¬ 
ciently large the usual Z' bounds are weakened by the necessity of producing it in association 
with b quarks. The dilepton resonance can also show up in decays of the B, in which case the 
final state would be a pair of dileptonic resonances and two b quarks, which can be paired up 
into two bii resonances. 

If (p is sufficiently heavy it can decay to Z'Z'. Or, if /) —)• Z'Z' is kinematically forbidden 
but the scalar mixing is sufficiently large, then p can decay to hh, WW, and ZZ if it is heavy 
enough. When B ^ pb dominates, BB production can therefore lead to events with as many 
as ten b quarks, with various sub-resonances among the 6-jets. Finally, if we incorporate DM 
into the theory the Z' and/or the p might decay invisibly, leading to events with 6-jets and 
MET. 
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Figure 5: BB production cross section at the LHC with ^/s = 14 TeV. 

Returning to our channel of primary focus, BB —?■ {bZ')(bZ') —>• b{bb)b{bb), the results of 
Section 4 are based on simulations of 45 separate parameter points covering a broad range of 
B and Z' masses. Before presenting those results, we describe our analysis technique. To aid 
in the discussion, we adopt three representative benchmark points. 

Benchmark 1 (Mb = 1 TeV, Mz' = 50 GeV) 

This light Z' benchmark is motivated by the secluded DM explanation of the GCE if, as 
discussed in Section 2.5, the DM mass is around 60 GeV. Larger values of Mz' are consistent 
with the GGE if the dark matter annihilates directly to bb. This benchmark requires jet- 
substructure techniques because the large mass difference between B and Z' means that the 
bb from the Z' decay will typically form a single massive jet. 

It is not difficult to find parameters consistent with Mb = 1 TeV, Mz' = 50 GeV, and 
Br{B —)• Z'b) ~ Br{Z' —?• 65) ~ 1. Eor example, start with g' = 1/20, corresponding to 
(<1>) = Mb/V2 and sr = A/\/2. Eor this value of g', B ^ cj)h is forbidden if the quartic 
coupling satisfies A$ > 1/2 (here we neglect scalar mixing), in which case Eigure 2 shows the 
B decays dominantly to Z'b (unless sr ~ 1). Eigure 3 shows that for k < 10“^A^, Z' will 
mainly decay to bb. If we incorporate Dirac fermion dark matter with = 60 GeV, the 
relic abundance requires g'Q^^ ~ 0.2 in the secluded DM scenario, or ~ 4. Then we need 
K < 3 X lO”"^ to satisfy direct detection constraints. 

Benchmark 2 (Mr = 1.5 TeV, Mz' = 750 GeV) 

This “medium mass” point can be discovered with high significance after 300 fb~^ of data, 
even with sizable systematic uncertainties, and will have hints after 30 fb“^ (see Eigure 10). 
An example set of model parameters for this point starts with ($) = 1500 GeV (corresponding 
to g' = 0.35 and sr = X). With this choice of parameters, > Mr is realized for A,^ > 1/4, 
in which case B —?■ Z'b typically dominates. Eor Z' —)• bb to dominate only requires k < 0.16A^. 
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Benchmark 3 {Mb = 2 TeV, Mz' = 1.5 TeV) 

Because of its small production cross section, this “high mass” point may require as much 
as 3000 fb“^ to be discovered. For an example set of parameters we can again start with 
(<!)) = 1500 GeV (corresponding to g' = ll\f2 and sr = 3A/4). To have > Mb we need 
A(^ > 4/9, and for Z' —)• bb to dominate we need k < 0.19A^. 


3.1 Simulation 

We implement the model in Feynrules [38]. Our signal simulations use MadGraph5_aMC@NLD 
[39] for parton-level event generation, PYTHIA_8.2 [40] for showering and hadronization, and 
DelphesS [41] for detector simulation. The dominant background comes from QCD mul¬ 
tijet production, followed by tt production. We simulate these background processes with 
PYTHIA_8.2 and DelphesS. Jets are clustered with Fast Jet [42] using the anti-kt algorithm 
[43] with R = 0.5. For Delphes settings we use the default “CMS” parameter card that comes 
with the distribution. This card sets the 6-tagging efficiencies for the high-py jets that will 
be important for our analysis at approximately 0.5 {\r]\ < 1.2) and 0.4 (1.2 < \r]\ < 2.5) for 
6-jets, 0.2 (jryj < 1.2) and 0.1 (1.2 < \ri\ < 2.5) for c-jets, and 10“^ for light jets. 

We use Hathor [44] to calculate vector-quark production cross sections at NNLO [45] with 
MSTW2008 NNLO parton distribution functions [46]. For the tt production cross section we 
take a = 954 pb, based on Ref. [45]. For the QCD background we adopt the LO cross 
section reported by Pythia. Pythia8 with default settings has been found to give reasonable 
agreement, at a level better than ~ 50%, with 7 TeV LHC data on multijet production 
[47, 48]. The difficulty in modeling the QCD background requires that it be estimated from 
data in an actual analysis. We discuss one approach to this estimation in Section 4. 

To reduce the statistical uncertainty associated with our QCD simulations, we bias the 
event generation to favor high-p^ events and record the event weights. We estimate the 
statistical uncertainties of our QCD Monte Carlo sample as 




(3.1) 


where the Wi are the individual event weights. This uncertainty is less than 10% for most of 
the signal windows we use to obtain the results of Section 4. 


3.2 Analysis 

Only jets with pT > 100 GeV and \r]\ < 2.5 are considered in our analysis. In the discussion 
that follows, “jet” refers to an object satisfying these criteria, and we calculate the scalar 
sum of jet pr’s. Hr, using only these jets. To be selected, an event must have at least four 
jets {rij > 4), three or more of which must be 6-tagged {rib ^ 3)- The probabilities to have 
various rib, among events with rij > 4 and Hr > 500 GeV, are shown in Figures 6 and 7 for 
the backgrounds and for the three benchmark points introduced above. 
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Figure 6: Probabilities to have 0, 1, 2, and 3 or more among background events with 

at least four jets and Ht > 500 GeV. For QCD events the probability to have at least 3 6-jets 
is 1.2 X 10“3. 


Signal Benchmark 1 Signal Benchmark 2 Signal Benchmark 3 



Figure 7: Probabilities to have 0, 1, 2, and 3 or more 6-jets, among background events with 
at least four jets and Ht > 500 GeV. The signal benchmarks are described in the text. 


Figure 7 shows a lower probability to satisfy the nj, > 3 requirement for Benchmark 1, 
because B decays produce highly boosted Z' particles for this parameter point, leading to Z' 
decays that typically produce a single jet. A more sophisticated analysis might attempt to 
keep track of the number of 6-tags associated with individual jets. Figures 6 and 7 also suggest 
that it may be advantageous to require more than three 6-jets, especially if one adopts a looser 
6-tag algorithm with a higher efficiency than we assume. For examples of how requiring a 
high number of b-tags (> 5) may be able to reduce backgrounds and allow discovery of certain 
signals, see Ref. [49]. We present results for an analysis based on nj, > 3 to be conservative, 
and we will see that with this analysis there is discovery potential for Mb = 2 TeV at the 
HL-LHC. 

For each selected event we apply three separate reconstruction strategies. These strategies 
differ in how many of the jets in the event are used in the reconstruction and in how Z' 
candidates are identified. Once Z' candidates are found the identification of B candidates 
proceeds identically for all three approaches. 

The four-jet reconstruction uses only the four hardest jets in the event. Among these 
four, two jets are identified as Z' candidates if their jet masses match to within 10% and 
both jets have r 2 /ri < 0.5, where tat is the V-subjettiness variable defined in Ref. [11] . This 
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Figure 8: Distribution of signal events satisfying rij > 4, n;, > 3, and Ht > 500 GeV, for the 
three benchmark signals. The rectangular selection windows are those of Table 1, optimized 
for 5 = 10%. 


approach is effective for Mb ^ Mz', in which case the Z' particles are produced with a large 
boost. The six-jet reconstruction uses the six hardest jets in the event. Among these six jets, 
two dijet pairs (comprising a total of four jets) are identihed as Z' candidates if the dijet 
masses match to within 10%. The five-jet reconstruction uses the hardest five jets and takes 
a composite approach. Among the hardest five jets, a single jet and a dijet pair are identified 
as Z' candidates if their masses match to within 10% and the single jet has r 2 /Ti < 0.5. 

Regardless of which reconstruction method is applied to a particular event, there remain 
two available jets after two Z’ candidates are identihed. These jets are paired with the 
Z' candidates in both possible ways. For each pairing, if the jet-Z' invariant masses are 
within 10% of each other, then the jet-Z' systems are identihed as B candidates, and the two 
{Mz',Mb) pairs are recorded. If Z' candidates cannot be used to hnd B candidates, then 
the Z' candidates are discarded along with their associated masses. 

A single event may yield numerous (M^/, Mb) pairs, produced by any and all of the three 
reconstruction methods. Once we establish a range of Mz' and Mb values as a useful window 
for a particular signal parameter point, we count an event as being in the window once and 
only once if any of its {Mz', Mb) pairs falls in that window. This single counting allows for 
a more straightforward statistical interpretation of results. 

Figure 8 shows the distribution of signal events in the M^/ — Mb plane for our three 
benchmarks. To make these plots we divide the Mz' — Mb plane into 10 GeV x 20 GeV 
pixels. A given event can count at most once in a given pixel but is allowed to be counted 
in multiple pixels. Similarly, Figure 9 shows the distribution of QGD and tt events in the 
Mz' — Mb plane. In the tt plot, we see a concentration of events near {Myi/,Mt) due to 
successful reconstruction of the W and top resonances. We get much larger counts in a bulk 
region whose position is set by the Ht, jet pT, and jet multiplicity requirements. These are 
combinatorially favored “mispairings” produced by the six-jet reconstruction. Mispairings 
also produce additional concentrations at Mg Mt, with Mz’ — My/ or with Mz' between 
My/ and Mt. Finally, the counts at very small values of reconstructed Mz' arise from the 
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four-jet reconstruction, where individual jets with similar small jet masses can constitute a 
pair of Z' candidates. 



Figure 9: Distribution of QCD (left) and tt (right) events satisfying nj > 4, n;, > 3, and 
Ht > 500 GeV. The pixels are 20 GeV x 40 GeV (left) and 10 GeV x 20 GeV (right) in 
size. The three signal selection windows from Figure 8 are also shown. 


The Z' jB search should be carried out in a way that covers as much of the Mz' — Mb 
plane as possible. In Section 4 we present the results of the following strategy: for a given 
Mz’-, Mb signal point to be tested, we impose a cnt Ht > 1.5Mb and construct an appropriate 
window in the Mz' — Mb plane. We set the boundaries of the window by a loose optimization 
of the quantity 

, ^ (3.2) 

y^S + B + (5B)^ 

where S and B are the expected numbers of signal and background events in the window, and 
where 6 represents the systematic uncertainty associated with the background in the window. 
Once a window is chosen, we quantify the expected significance of the signal using 


S 

+ (6B)^ 


(3.3) 


for B > 50. For smaller B, we take the background-model probability of observing n counts 
as 

P(n) = J dXf(X\B,SB)g(n\X), (3.4) 

where g(n\X) is the Poisson distribution with mean A and f(x\^, a) is the normal distribution 
with mean ^ and standard deviation a, and we quantify the significance as 


2 ^ P(n) 

n>S+B 


erf ^ 


1 


(3.5) 























Benchmark 1 
[50 GeV, 1000 GeV] 

Benchmark 2 
[750 GeV, 1500 GeV] 

Benchmark 3 
[1500 GeV, 2000 GeV] 

Bottom-left corner 
Top-right corner 

(30, 840); (40, 860) 
(70, 1120); (60, 1100) 

(640, 1280); (660, 1360) 
(780, 1560); (780, 1520) 

(1360, 1840); (1360, 1840) 
(1500, 2080); (1500, 2080) 


Table 1: Benchmark signal windows optimized for L = 300 fb with all units in GeV. 
The first entry is for 0% systematics (h = 0), and the second entry is for 10% systematics 
(h = 10%). 


We will present results for <5 = 0 and 5 = 10%. In the following section, we argue that 
estimating background from data at a 10% level or better is a realistic goal for this analysis. 


4 Results 


We have studied the discovery prospects for 45 signal parameter points in all. Tables 1 and 
and 2 provide detailed results for our three benchmark points. Table 1 gives the Mz' — Mb 
selection windows used for each benchmark, optimized for an integrated luminosity of £ = 300 
fb“^, and for either h = 0 and 6 = 10%. The windows for 6 = 10% are shown in Figures 8 
and 9. Table 2 shows the numbers of events that pass the various cuts in our analysis, for 
background and for the three signal benchmarks. 

In an actual experimental analysis it will be important to estimate the QCD background 
from data. The background in a given window can be estimated using events with fewer 
6-tagged jets. For the <5 = 10% selection windows of Table 1, Table 3 compares the numbers 
of background events that pass the full analysis with the estimate 


E 

nj>4 


. rij jets and > 3\ 

74.,- ipf.s p.nH n.i. = 0 / 


X in window, with nj jets and rib = 0) 


(4.1) 


In the first factor, the events must pass the Ht cut (which differs for the different benchmarks, 
as the Ht cut is set to be Ht > 1.5Mb), but the events are not required to yield Z' or B 
candidates. In the second factor, the events must pass the full analysis, with at least one pair 
of Z' and B candidates with masses in the window, except that the usual requirement re;, > 3 
is replaced with re^ = 0. 

Instead of using re^ = 0 events for the estimate, one could instead use rib < 3, re^ = 1, 
or rib = 2 events, which might be more accurate. However, Table 3 shows that using re^ = 0 
events works rather well for the benchmark windows, and the signal contamination of the 
background in the re^ = 0 samples is less than 1% for all three windows. 

For most of the signal points we investigated, the accuracy of the estimate using rib = 0 
events is comparable to the level of agreement shown in Table 3. Exceptions include several 
of the points with M^/ = 100 GeV, where the tt background makes up a larger component of 
the background then for other points, due to the presence of W's. However, these points are 
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Cut 

4 jets, Ht > 500 GeV 
Ht > 1500 GeV, Ub > 3 
Ht > 2250 GeV, nb>3 
Ht > 3000 GeV, nb>3 

Mass pair within 10% 


Sig. 1 analysis 0% (10%) 
Sig. 2 analysis 0% (10%) 
Sig. 3 analysis 0% (10%) 


QGD 
8.9 X 10® 
47200 
5550 
834 

/6850\ 
1030 

V158/ 

41.8 (18.3) 
130 (72.5) 
10.5 (10.5) 


tt 

6.9 X 10® 
5400 
643 
98.4 
/1080\ 
132 

4.34 (1.67) 
15.0 (8.35) 
1.09 (1.09) 


Benchmark 1 
14900 
3740 
1160 
203 
/ 644\ 
260 

V54V 

276 (256) 


Benchmark 2 Benchmark 3 




109 (81.1) 


5.51 (5.51) 


Table 2: Cut table for C = 300 fb“^. In the fifth row the three entries are for cuts on Ht 
of 1500, 2250, 3000 GeV, respectively. The final results, shown in the bottom three rows for 
(5 = 0 and 6 = 10%, have Ht > and require events to land in the appropriate {Mz’, Mb) 
window, defined in Table 1. 



Benchmark Window 1 

Benchmark Window 2 

Benchmark Window 3 

full analysis, with Ub >3 

19.9 ±1.1 

80.9 ± 1.7 

11.6 ±0.3 

rib = 0 

rib >3 estimate 
S/B in window 

21.7 ±0.4 ±0.3 

8.3 X 10“3 

78.5 ± 0.5 ± 1.9 

3.5 X 10“^ 

12.0 ±0.1 ±0.7 

1.7 X 10“^ 

rib = 1 

rib >3 estimate 
S/B in window 

21.9 ±0.4 ±0.6 

0.15 

79.6 ±0.6 ±2.4 

6.3 X 10“3 

12.1 ±0.1 ±0.9 

2.9 X 10"3 

rib = 2 

rib >3 estimate 
S/B in window 

21.2 ±0.6 ±1.5 

1.4 

79.7 ±0.9 ±4.0 

6.4 X 10“2 

12.0 ± 0.2 ± 1.6 

2.9 X 10“2 


Table 3: For £ = 300 fb“^, actual background counts (top row) and the associated estimates 
using events with zero, one, or two 6-tagged jets. The actual counts come with Monte Carlo 
uncertainties, and the estimates come with Monte Carlo uncertainties followed by statistical 
uncertainties associated with the estimation method. Also shown are signal-to-background 
ratios for each window and Ub requirement. 


heavily signal-dominated, i.e. they have a large S/y/~B. If the background estimation is not 
quite as good as we assume, the discovery potential changes very little. Furthermore, other 
handles for estimating the background will be at experimentalists’ disposal, including events 
with reconstructed Mz' and/or Mb values outside the window, or perhaps events for which 
the mass-matching that identifies Z' and/or B candidates fails at 10% but satisfies some less 
stringent requirement. 

Figure 10 shows the projected discovery potential in the Mz' — Mb plane for £ = 30 
fb“^, 300 fb“^, and 3000 fb“^. Discovery at the 5a level is possible for a broad range of Mz', 
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Figure 10: Significance contonrs for £ = 30, 300, and 3000 fb and for 5 = 0 and 10%. The 
black points indicate the signal parameter points from which the contonr plot was generated. 

with Mb < 1250 GeV for 30 fb“^, with Mb < 1600 GeV for 300 fb“^, and with Mb < 2000 
GeV for 3000 fb-\ 

5 Conclusions 

The hnnt for new colored fermions is an integral part of the broad search strategy employed 
at the LHC. To date almost all searches for new vector-like partners of the top or bottom 
quarks have been in final states containing SM bosons {W, Z, or h). We have pointed out 
that, by virtue of being vector-like, it is straightforward for the heavy qnarks to be charged 
under additional gauge groups, and that these couplings may dominate their decays. We 
have focussed on the simple case of a new U{iy gronp which a vector-like B quark is charged 
under. We have described a simple realisation of this scenario, based around the concept of the 
“Effective Z'”. We have outlined the wide range of new phenomena and interesting search 
channels that exist in this class of simple models, which contain only three new particles. 
If the kinetic mixing between U{iy and hypercharge is small the new channels all involve 
mnltiple h qnarks. We demonstrated that there is a broad region of parameter space in these 
models where the new decay B —)• hZ' —)• b{bb) dominates. 

We have presented a search method that can simultaneously observe the new quark and 
the new gauge boson in final states containing np to six b quarks, by carrying out a two- 
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dimensional mass reconstruction of events. The large QCD and smaller tt backgrounds can 
be effectively reduced by requiring pairs of resonances whose masses are close, which in turn 
contain sub-resonances whose masses reconstruct to be the same. Although there are many 
b quarks in the final state we take a conservative approach and require only three 6-tags. A 
better understanding of 6-tagging efficiencies may allow this requirement to be strengthened, 
leading to a further suppression of background. The kinematics of the process are sensitive 
to the mass splitting between B and Z' and we account for this be varying our reconstruction 
technique with the number of final state jets and employing the techniques of A^-subjettiness 
to uncover merged jets from the Z' decay. We find that discovery at the 5a level is possible 
for a broad range of Mz', with Mb < 1250 GeV for 30 fb“^, with Mb < 1600 GeV for 300 
fb“^, and with Mb < 2000 GeV for 3000 fb“^. 

It is intriguing that the recently observed Galactic center excess can be explained by 
weak scale DM annihilating into 6 quarks. If this takes place through a new mediator one 
might expect new 6-quark partners which may themselves decay into the mediator. We have 
provided one such example of this and have shown that the LHC has the capability to test 
this DM scenario over much of its parameter space. 

Finally, the technique we describe is not unique to the model we analyse and will be 
widely applicable to many models where a new particle is pair produced and decays to a 
lighter new state, finally decaying to SM particles. For instance, the approach we advocate 
has an obvious extension to vector-like top quarks, T —)• tZ' —)• t{bh)/{tt). It would also 
enhance RPV gluino searches [50, 51] in the case where the squarks are lighter than the 
gluinos. 

Note Added 

While this work was in the final stages of completion CMS released details of a search for T in 
the exotic mode T —)• bW' with W decaying leptonically [52]. The CMS analysis also searches 
simultaneously for two new particles and carries out a two-dimensional mass reconstruction 
of events, but the final state and particle content are different from what we consider. 
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